Enhanced Stability of Bound Pairs at Nonzero Lattice Momenta 
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A two-body problem on the square lattice is analyzed. The interaction potential consists of strong 
on-site repulsion and nearest-neighbor attraction. The exact pairing conditions are derived for s-, 
p-, and d-symmetric bound states. The pairing conditions are strong functions of the total pair 
momentum K. It is found that the stability of pairs increases with K. At weak attraction, the pairs 
do not form at the F-point but stabilize at lattice momenta close to the Brillouin zone boundary. The 
phase boundaries in the momentum space, which separate stable and unstable pairs are calculated. 
It is found that the pairs are formed easier along the (n, 0) direction than along the (n, n) direction. 
This might lead to the appearance of "hot pairing spots" on the K x and K v axes. 

PACS numbers: 03.65.Ge, 71.10.Li 
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I. INTRODUCTION 

The short coherence length observed in high-tempera- 
ture superconductors has stimulated theoretical research 
on real-space pairing on the lattice. Strong electron- 
phonon and electron-electron interactions coupled with 
week screening and reduced dimensionality lead to for- 
mation of preformed bound pairs that Bose-condense and 
form a superfluid stated. The two interacting particles 
on a lattice constitute a class of exactly solvable quan- 
tum mechanical problems. A number of physical charac- 
teristics can be obtained in closed analytical forms: the 
pairing threshold, binding energy, wave function, effec- 
tive radius, and others. It is essential that the reduc- 
tion of the two-body Schrodinger equation to a one-body 
Schrodinger equation occurs differently than in a contin- 
uous space. The underlying lattice introduces a preferred 
reference frame, which results in non-trivial dependencies 
of the physical quantities on the total lattice momentum 
of the pair K. 

In the present work we show that bound states on a 
lattice become more stable at large K. This property is 
opposite to the Cooper pairs in a BCS superconductor. 
The effect is best seen through the K-dependence of the 
threshold value of the attractive potential. In two di- 
mensions, the threshold is always zero if the potential is 
purely attractive. Therefore, it is essential to consider 
mixed repulsive-attractive potentials. (Such potentials 
are also more realistic.) A natural choice is the Hub- 
bard on-site repulsion U combined with an attraction V 
acting on one or more layers of nearest neighbors. The 
latter mimics the overscreened Coulomb interaction re- 
sulted from electron-phonon interactions with distant lat- 
tice ions&i. 

General properties of the two-body states on a lattice 
were reviewed by Mattis*. Lin gave the first solution 
for the s-symmetric singlet ground state in the U — V 
model 6 . In fact, that solution was presented as the low- 
density limit of the t — J model. Pctukhov, Galan, and 
Verges extended the solution to the p-symmetric triplet 
and d-symmetric singlet pairs, again in the framework 
of the t — J model. Kornilovitch found pairing thresh- 



olds for a family of U — V models, in which attrac- 
tion extended beyond the first nearest neighbors^. Basu, 
Gooding, and Leung studied pairing in a U — V model 
with anisotropic hopping*!. Two-body problems were also 
solved for several models with multi-band single-particle 

Previous studies were mostly confined to the T-point of 
the pair Brillouin zone. (Reference does contain some 
results for the diagonal and the boundary of the Brillouin 
zone.) In this work we focus on the K-dependence of the 
pairs properties. We formulate and solve the pair stabil- 
ity problem in its general form. We define the pairing 
surface inside the Brillouin zone, which separates the re- 
gions of pair stability from the regions of pair instability. 
We also show that, in general, the shape of the pairing 
surface is different from the shape of the single-particle 
Fermi surface. The difference leads to the appearance of 
pairing "hot spots" , which are regions in the Brillouin 
zone where real space pairing takes place while being en- 
ergetically unfavorable elsewhere. 



II. THE TWO-PARTICLE SCHRODINGER 
EQUATION 

In this paper we consider the simplest U — V model on 
the square lattice, the one that includes attraction be- 
tween the first nearest neighbors only. For more complex 
attractive potentials, see Ref. The model Hamilto- 
nian reads 



H 



nj 



V ^2 n m n r 



(1) 



Here m = n + 1 denotes a lattice site that is a near- 
est neighbor to n, n na — c^Cnj, and n n = n n f + n n |. 
The rest of notation is standard. We allow one spin- 
up fermion and one spin-down fermion in the system. 
It is convenient to not to fix the permutation symme- 
try of the coordinate wave function from outset. In this 
way, the singlet and triplet solutions will be obtained 
simultaneously. The system is fully described by the 
tight-binding amplitude ^(ni, 112) or, equivalently, by its 
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Fourier transform $(ki,k 2 ). Fourier transformation of 
the Schrodinger equation with Hamiltonian yields an 
equation for <!>: 

[E - £ (k 1 )-£(k 2 )]$(k 1 ,k 2 ) = f/^$(q,k 1 +k 2 -q) 

q 

■ ^E e " kll E $ ^ k i+ k 2- c i) e4q1 ' ( 2 ) 



where 1 = (±1, 0), (0, ±1) are the four nearest-neighbor 
vectors of the square lattice, E is the combined en- 
ergy of the particles, and e(k) = — tJ2\ ex P(~ *kl) = 
— 2t(cosk x + cosfcj,) is the free single-particle spectrum. 
The integral equation @ is solved by introducing five 
functions A(K) that depend only on the total lattice mo- 
mentum K = ki + k 2 : 



A (0 ,o)(K) = 2^ *(q, ki + k 2 q), 
q 

AiCKjs^^q.ki+ka-qJe*'* 1 . 
q 

The wave function is then expressed as follows: 



*(ki,k 2 ) = 



t/A (0 , 0) (K)-^EiAi(K) t 
E - e(k x ) - e(k 2 ) 



-ikil 



(3) 
(4) 

(5) 



Substituting this solution back into the definitions @ 
and Q one obtains a system of linear equations for 
A( ,o), A (+i,o), A (-i,o), A (o,+i), and A (0 The de- 
terminant of that system defines energy E as a function 
of K: 
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where 



L P =L P (E,K) = J2 



E - e(q) - e(K - q) 



(6) 



(7) 



The eigenvector of the above system, upon substitution 
in Eq. JSJ , defines the pair wave function up to a normal- 
ization constant. 

The derived formulae illustrate the general structure 
of two-body lattice solutions. The spectrum equation 
is a determinant, with the dimension equal to the to- 
tal number of sites within the interaction range. In the 
present model this number is five, one for the Hubbard 
term and four for the nearest-neighbor attraction. An 
eigenvector of the corresponding matrix determines the 
two-particle wave function. The matrix elements are mo- 
mentum integrals of the ratio of two energy polynomials 
P n 2_ 1 (E)/P n 2 (E), where n is the number of one-particle 
energy bands. On the square lattice n = 1, hence the 



simple form of the integrals L p . The copper-oxygen plane 
is a more complex case with n = 3. The corresponding 
two-body lattice problem was solved in Ref. [Til ]. 

Analytical investigation of the spectrum and wave 
function may be difficult, especially in cases when matrix 
elements are not calculable in closed forms. The pairing 
threshold is easier to find because the determinant needs 
to be analyzed only at one special energy point. At the 
same time, the threshold is a convenient and sensitive in- 
dicator of pair stability. The binding condition of model 
Q will be the prime focus of the rest of this paper. 

The spectrum equation (JSJ) depends on three parame- 
ters: the interaction potentials U and V, and the total 
lattice momentum K of the pair. A binding condition is 
defined as a functional relation between U, V, and K, 
when the energy of a bound state becomes equal to the 
lowest energy of two free particles with the same com- 
bined momentum K. At fixed K, the minimum energy is 
realized when both particles have equal momenta K/2. 
This is because the kinetic energy of relative motion is 
zero. This statement can also be derived rigorously by 
minimizing a sum of two single-particle energies. The 
minimum energy is given by 

Eb(K)=2e(y) =-4icos^-4icos^. (8) 

Thus real-space bound states are formed by particles 
moving parallel to each other. This should be compared 
with the Cooper pairs, which are formed by quasiparti- 
cles with opposite momenta. 

Imagine a bound state with energy E(K) just below 
£b(K), and let E(K) — > E (K). The spectrum equation 
then defines a relationship between U, V, and K, which 
is the pairing condition of interest. Before attempting 
general solution let us consider the T-point. At K = 
(0, 0) the solution acquires additional symmetry, which 
simplifies analysis of the spectrum equation. It follows 
from definitions JJJ that L qx = L_ q% — L q = L— q , 

" L n 



L-a. 



2q. T 
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L-q x -q v 



2<?, y — L-2q y , and L qx+qy — ±^ q ^- qy 
Next, introduced a new basis: 



Ho,o)j 



A 



[A(i,o) + A(_i i0 ) + A( ,i) + A( 0l _i)], 

A (o,-i)L 



h- - t [A(i,o) - A(-i,o) + A( ,i 



A p _ 



Ho, 



i)]. 



IA(i,o) - A(-i,o) - A (0 !) 
IA(i,o) + A(-i,o) - A (CU ) - A( 0i _i)]. (9) 



The new basis functions are eigen-functions of the sym- 
metry operators of the square lattice. States of differ- 
ent symmetries separate and matrix (0 becomes block- 
diagonal. We consider the blocks separately. 

s- symmetry. The 2x2 block that mixes Ao and A s 
has the form 



Li 



on 
L' 



L'r 



L' 



A 
A, 



= 0, 



(10) 
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where 



L 'oo ~ E 



L' 



E 



L 'ss = E 



E - 2e(q) ' 

cosg a 
E - 2e(q) ' 

2 cos (fc (cos (fc + cos q y ) 
E - 2e(q) • 



(11) 



and Ln, 



4L' n . Substitution E = En 



-8t leads to 



logarithmic divergence of all integrals L' . In order to cir- 
cumvent this difficulty, add and subtract Lq , 4Lq , and 
4L' 0Q from L' sQ , L' 0s , and L^, 



ences converge in the limit E 



respectively. The differ- 

— > —8t'. L.n — Lnn = at, 



J 0s 



4 ^oo — 2t' anc ^ - 4 ^oo — 2t' 



J s0 -^00 — 8t' 

Upon expansion 



of the 2x2 determinant the (L' 00 ) 2 terms cancel while 
the rest yields 



(2tU - 8tV - UV) ■ L' m 



V + 2t 



UV 



0. (12) 



The equation is satisfied when the coefficient at the di- 
vergent Lq is zero. This results in the binding condition 
for the s-symmetric paii-Si 



V > V, = 



2Ut 
U + 8t 



(13) 



In the weak coupling limit U, V <C t, the threshold re- 
duces to V s = jU. The factor 4 here is the number of 
nearest neighbors: the on-site repulsion needs to be 4 
times stronger to balance the attraction acting on 4 sites 
at once. In the limit of infinite Hubbard repulsion, the 
attraction threshold V S (U — > oo) = 2t is finite. Limited 
influence of an infinite repulsive potential is explained by 
vanishing of the wave function at ri = rj. Without the 
repulsion, the attractive threshold goes to zero, as ex- 
pected in two dimensions. Thus the Hubbard term raises 
V s from zero to 2t. 

p- symmetry. Two lxl blocks describe two-particle 
states with p orbital symmetry. Their energy is deter- 
mined by the equation 



E 



o • 2 

z sin q x 
E - 2e(q) 



1 

V 



0. 



(14) 



Unlike the s-symmetry case the integral converges at E = 
— 8t. Integration results in the binding condition for p- 
symmetric pairsi 



2tt 

V>V v = —2 



t = 5.50 1. 



(15) 



Notice that neither the spectrum nor the binding condi- 
tion depends on the Hubbard potential U. 

d-symmetry. The last lxl block describes two-particle 
states with d orbital symmetry. The spectrum equation 
is 



E 



2 cos q x (cos q x — cos q y ) 
E - 2e(q) 



(16) 



s and d singlets and two p triplets 



one s singlet and two p triplets 



one s singlet bound state 




FIG. 1: Phase diagram of the two-particle model fl) at K = 
(0,0). No bound states exist below the lowest solid curve. 
The phase boundaries are described by Eqs. \VA) . I|15|l . and 
117H . The dashed line indicates the limit value of the s-state 
threshold V s = It for an infinitely strong Hubbard potential. 



Again, the integral converges in the limit E 
leads to the binding condition 7 : 



V>V d = 



2tt 



t= 7.32 1. 



-8t, which 



(17) 



As with p-symmetric pairs, the binding condition does 
not depend on U. This is because the wave functions of 
the p and d states are identically zero at ri — T2- The 
contact potential is not felt even at small V. 

Figure^shows the "phase diagram" of the two-particle 
system at the T-point. 



III. 



THE GENERAL BINDING CONDITION 



Let us return to the full spectrum equation (|SJ| and find 
the binding condition at an arbitrary pair momentum K. 
First, using the explicit form of e(q) and shifting q by 
K/2 the integrals L p are transformed as follows: 



L 



L 
L qx 

q*+q y 



E 



i 



E + a cos q x + b cos q x 



= M, 



on- 



= p.' — 



E 

q 

*E 



cos q x 



E + a cos q x + b cos q x 



e l ~M 



10 j 



JK 



q 

-K^+Ky 



cos 2q x 



E + a cos q x + b cos q x 



Al- 



io- 



E 



cos q x cos q y 



E + a cos q x + b cos q x 



■ K X +Ky 

e *^s— Mil, 



and so on, where 



a = Atcos(K x /2) > 
b = 4tcos(K y /2) > 0. 



(18) 



(19) 
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The subscripts of M nm refer to the frequencies of cos nq x 
and cosmq y in the numerator of the integrand. Upon 
substitution E — > E (K) all M nm diverge logarithmi- 
cally. Similarly to the T-point case, a divergent Moo is 
subtracted and added to each M nm . The differences con- 
verge. A straightforward calculation yields: 



Mi,, 
Moi 
M 20 

M 02 
Mil 



M10 - M 00 = — arcsin J — — , 
ira V a + 

M i - M 00 = —7 arcsin - 

no V a + 



4 / a + & 



■ arcsm 



ira \ a 



a + b 



4 a + b b 

— — - — arcsin \ - 

nb \ b V a + b 



M xl - Moo = 




On the next step, the determinant JBJ is expanded result- 
ing in a lengthy expression, which is in general a fifth- 
order polynomial in Moo- However, the coefficients at 
the second through fifth powers cancel identically. The 
spectrum equation reduces to A ■ Afoo + B = 0, where 
A and B are some complicated functions of M nm . [A is 
given below as the l.h.s. ofEq. 1)21(1.] Cancellation of the 
high-order powers of logarithmically divergent terms was 
observed in other two-dimensional modelE&iLiS. It seems 
to be a general property of the pairing problem in two 
dimensions, although no rigorous proof is given here. In 
the present model, the divergence of M o implies A = 0. 
After lengthy transformations and factorizations the last 
condition can be written as follows: 



(20) 



{1-VM 20 }- {1-VMoa}- 

{UV 2 [(Mao - 4M 10 )(M 02 - 4M 01 ) - 4(M n - M 10 - M 01 ) 2 ] 

+ UV [(Mao - 4Mi ) + (M 02 - 4M i)] + U 2 [8M U - 2(M 20 + M 02 )] + U - W) = 0. (21) 



This is the general binding condition sought. It repre- 
sents the main analytical result of the paper. The first 
two factors correspond to the two p-symmetric triplet 
pairs. This is best seen at the diagonal of the Bril- 
louin zone K x = K y , which implies a — b. This yields 
M 02 = M 20 = (tt- 2)/(2ntcos(K x /2)), and the binding 
conditions become 



V > V P (K X = K y ) 



2tt 



tcos (K x /2) 



7T-2 

5.50 1 cos (K x /2). (22) 



At K x = 0, the last equation reduces to Eq. 1)15(1. Hence, 
it is identified with the p-states. 

Equation (|22|l is a particular example of the general 
effect: the enhanced stability of bound states at large 
lattice momenta. Indeed, the pairing threshold decreases 
with K x and vanishes completely in the corner of the 
Brillouin zone. At any V < 5.50 1, there exists a pair- 
ing surface in the Brillouin zone, which separates stable 
and unstable bound states with p orbital symmetry. The 
concrete form of the pairing surface can be found numer- 
ically. 

Next, consider the last factor in Eq. I|21|l . which de- 
scribes the s- and d-symmetric pairs. The binding condi- 
tion does not factorize at arbitrary K. However, factor- 
ization occurs at the diagonal of the Brillouin zone where 



M 02 = M 20 and M i = M w = [8tcos(K x /2)}- 1 

s : UV(M 20 + 2M U - 8M10) + U - 4V = 0, 
d : V(M 20 - 2M n ) + 1 = 0. (23) 

Substitution of the integrals from Eq. ()20|l results in the 
following binding conditions for the s and d states: 



V > V S {K X =K y ) 



V > V d (K x = Ky) 



2Utcos (K x /2) 
U + 8tcos(K x /2Y 



2n 



icos (K x /2) 



4 - 7T 

7.32*008(^/2) 



(24) 



(25) 



Comparing the binding conditions on the Brillouin zone 
diagonal with the T-point, one observes that in both cases 
the hopping integral t is replaced with tcos (K x /2). The 
s binding condition approaches 2t cos (K x /2) in the limit 
of infinite U. The d binding condition does not depend 
on U at all. Again, the stability of bound states increases 
with the pair momentum. 

At an arbitrary pair momentum, the pairing surface 
can be found numerically by solving the transcendental 
equation 1)21(1. Figure |2 shows the pairing surfaces for 
V = 1.5 1 and U = 50 i. If < V < 2Ut/(U + St), there 
are four pairing surfaces corresponding to the s, two p, 
and d pairs. At V = 2Ut/(U + 8t), the s pairing sur- 
face shrinks to a point indicating that the s-symmetrical 
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FIG. 2: Pairing surfaces for U = 50 i and V = 1.5 1. Bound 
pairs are stable outside the respective surfaces. At these pa- 
rameters both p and d surfaces terminate at the Brillouin 
zone boundaries: p at (±0.46 it, ±7t) and (±7r, ±0.46 7r); d at 
(±0.41 7r, ±7r) and (±7r, ±0.41 7r). The p and d surfaces actu- 
ally cross. 



bound pair is the ground state of the system at any mo- 
mentum K. Similarly, there are two p and one d pairing 
surfaces at 2Ut/{U + 8t) < V < 5.50 1, only one d pairing 
surface at 5.50 1 < V < 7.32 1, and no surfaces at all at 

V > 7.32*. 

Let us discuss the evolution of pairing surfaces with 

V in more detail. Figure [21 shows the evolution of the s 
surface at U = 50 1. At small V, the pairing surface ap- 
proaches the boundaries of the Brillouin zone. For most 
K, the strong repulsion overcomes the weak attraction 
and the ground state consists of two unbound particles. 
But in the very vicinity of the Brillouin zone boundary 
the scattering phase space shrinks, allowing formation of 
a bound state with a binding energy « V. The pair sta- 
bility region exists for any finite V, however small. The 
pairing surface never crosses the boundary of the Bril- 
louin zone. (The latter fact was noticed in Ref. Q.) As 

V increases the stability region expands. The pairing 
surface shrinks while assuming a circular shape. Finally, 
at V = 2Ut/(U + St) it collapses to the T-point. 

Figure 01 shows the ^-evolution of one of the p-states 
pairing surfaces. At small V the surface terminates at 
some points along the K y = ±n boundaries. The K x 
coordinate of the end points is given by the condition 
2tcos(K x /2) = V. There are no end points along the 
K, x = ±7r boundaries. At V = 2t, the two pairs of 
end points merge at (0, ±7r) resulting in an elliptically 
shaped critical surface. As V increases further, the crit- 
ical surface shrinks to the origin while retaining the el- 
li ptic shape. The r atio of the two axes of the ellipse is 
\J (3-7T — 8)/(4 — 7r) = 1.29. The pairing surface for the 



FIG. 3: The pairing surfaces of the s-symmetrical bound 
state for U = 50 1 and several V , obtained by solving Eq. 11' U . 
Going from outside inward, the plots correspond to V = 0.2 1, 
0.5 1, 1.0 1, 1.3 1, 1.5 1, 1.6 1, and 1.7 1. The surface shrinks to 
the origin at V = 2Ut/(U + 8t) = 1.724 1. The pairing surface 
never crosses the boundaries of the Brillouin zone. 
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FIG. 4: The pairing surface of one of the two p-symmetric 
bound states for U — 50 1. Going from outside inward, 
the lines correspond to V = 1.0 1, 2.0 1, 3.0 1, 4.0 1, 5.0 1, 
and 5.4 1. The four end points of the V = 1.0 1 plot are 
indicated by semicircles. When V — > 5.50 1, the surface 
shrinks to the origin but remains elliptic with the eccentricity 
e = 2y (tt — 3)/(37r — 8) = 0.63. The second family of the p 
pairing surfaces is obtained by rotating the plots by ±7r/2. 
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FIG. 5: The pairing surface compared with a single particle 
Fermi surface. The thin solid line is the Fermi surface defined 
by the condition e(kp) = —2.1t. The thick solid line is the 
same Fermi surface multiplied by 2. The dashed line is the s- 
symmetric pairing surface for U — 50 1 and V — 0.8 1. Notice 
the "binding hot spots" along the x and y axes where the 
double Fermi surface extends beyond the pairing surface. 



second p-symmetric bound state is obtained by rotating 
the plots of Fig.Hby ±n/2. 

The d-state pairing surface is symmetric under the 
K x <-* K y exchange, as evident from Fig. [21 There 
are four pairs of end points with the non-trivial mo- 
mentum coordinate determined by the condition 2t(U — 
2V) cos {K x . v /2) = UV. At V = 2Ut/{U + At), the end 
points merge at (0, ±7r) and (±7r,0). At larger V the 
pairing surface assumes the circular shape similar to the 
s-surface shown in Fig. [3] The d-surface shrinks to the 
origin at V = 7.32 1. 



IV. DISCUSSION 

In continuous space the internal (binding) energy of a 
composite non-relativistic particle is independent of its 
momentum, and the effective mass is independent of the 
internal energy. In contrast, on a lattice the total mass 
increases with the internal binding energy 5 . In this pa- 
per, we established another unusual property of lattice 
pairs: the stability of a bound state increases with its 
momentum. Although the total pair energy grows with 
the momentum, the corresponding energy of two free par- 
ticles grows faster. As a result, at the Brillouin zone 
boundary the energy of the bound state always lies below 
the continuum of single-particle states. At the boundary, 
the periodicity of the pair wave function is commensurate 
with the lattice. The particles could choose to coherently 
"see" only one component of the potential. At the same 



time, the kinetic energy of the relative motion is sup- 
pressed because of shrinking of the phase space available 
for mutual scattering. Thus as long as the interaction 
potential has at least one attractive region the particles 
will bind. 

In two dimensions the effect is somewhat hidden by 
the fact that an arbitrary weak attractive potential binds 
particles even at the T-point. As K increases no quali- 
tative changes happen and the pair remains stable. The 
situation is different when the potential has strong repul- 
sive pieces. At the T-point the repulsion overcomes the 
attraction resulting in a non-bound ground state. But at 
the Brilloiun zone boundary the ground state is always 
bound. By continuity, a pairing surface must exist that 
separates stable from unstable bound states. The shape 
and size of the pairing surface is a sensitive characteristic 
of the potential. Examples of such surfaces were given in 
the previous section. 

Imagine now a finite density of fermions. As filling 
increases, particles with large momenta become avail- 
able. At some filling, the double Fermi surface crosses 
the pairing surface. (The single-particle Fermi surface 
should be multiplied by a factor of 2 for adequate compar- 
ison.) In general the two surfaces have different shapes. 
Therefore crossing begins at some "hot" segments of the 
Fermi surface. The phenomenon is illustrated in Fig- 
ure |SJ The pairing surface for U = 50t and V = 0.8 1 
is compared with the double Fermi surface correspond- 
ing to a Fermi energy —2.1 1 (filling = 0.35 particles per 
unit cell). Notice how the pairing surface is elongated 
along the diagonals of the Brillouin zone while the Fermi 
surface is elongated along the K x and K y axes. The 
overlay of the surfaces creates four regions near (±7r, 0), 
(0, ±7r) where pairing correlations are enhanced. Un- 
der such conditions a particle with a momentum close 
to (±7r/2,0) or (0, ±7r/2) can pick up a second particle 
with the same momentum and form a bound state. Such 
a pairing mechanism is in contrast with the Cooper pair- 
ing. Obviously, the presence of filled states modifies the 
low-energy scattering dynamics of the pair. Analysis of 
the many-body model Q near a paring instability is a 
difficult problem, which warrants a separate investiga- 
tion. An interesting question is whether the ll U -negative 
V" model possesses a state where bound pairs with non- 
zero momenta "float" atop a fermion sea. Apart from 
academic interest, such a state may be of interest for the 
pseudogap phenomenon in high-temperature supercon- 
ductors. (For the latter see, e.g., Ref. 0.) One possible 
approach to the problem would be the scattering matrix 
approximation, which has been successfully applied to 
the negative Hubbard modeli^i 5 -. 

In conclusion, we have considered a two-particle lattice 
problem with on-site repulsion and nearest-neighbor at- 
traction. The stability of the bound states increases with 
the pair's total momentum, which is a consequence of the 
lattice discreteness. The effect is conveniently described 
in terms of the pairing surface that separates stable from 
unstable pairs. Exact pairing surfaces for s-, p-, and d- 
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symmetric bound states have been found and analyzed. possible, caused by the difference in shapes between the 
At finite fillings, the appearance of "hot pairing spots" is pairing and Fermi surfaces. 
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where 



L p = L p (E,K)^J2 



E - e{q) - e(K - q)' 



(A2) 



and e{k) = —2tcosk. The matrix elements are readily 
calculated. Expansion of the determinant then yields two 
solutions. 

Positive spatial parity (singlet bound state). The sin- 
glet energy E s < is determined from the equation 



yjEl - 16 i 2 cos 2 {K/2) + Uj ■ 

■ (yE 2 s -l6t 2 cos 2 {K/2) - E t 



-2V{U - E s ) = 0. (A3) 



The bound state stabilizes at 



V > 



2Utcos(K/2) 
U + Atcos{K/2)' 



(A4) 



APPENDIX A: BOUND STATES IN ONE 
DIMENSION 



Negative spatial parity (triplet bound state). The 
triplet energy is given by 



In this appendix we present, for reference, the two- 
body solution of model JQ) is one dimension. The spec- 
trum equation is a 3 x 3 determinant 



4t 2 

E t = -V-— cos 2 {K/2). 



(A5) 



u 



-2q 



0. 



The triplet is stable when V > 2tcos{K/2). Notice that 
(Al) in the U — > oo limit, E s = E t . 



